Emergence of grouping in multi-resource minority game dynamics 
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The Minority Game (MG) has become a paradigm to probe complex social and economical phenomena where 
adaptive agents compete for a limited resource, and it finds applications in statistical and nonlinear physics as 
well. In the traditional MG model, agents are assumed to have access to global information about the past 
history of the underlying system, and they react by choosing one of the two available options associated with 
a single resource. Complex systems arising in a modem society, however, can possess many resources so that 
the number of available strategies/resources can be multiple. We propose a class of models to investigate MG 
dynamics with multiple strategies. In particular, in such a system, at any time an agent can either choose a 
minority strategy (say with probability p) based on available local information or simply choose a strategy 
randomly (with probability 1 — p). The parameter p thus defines the minority-preference probability, which is 
key to the dynamics of the underlying system. A striking finding is the emergence of strategy-grouping states 
where a particular number of agents choose a particular subset of strategies. We develop an analytic theory 
based on the mean-field framework to understand the "bifurcation" to the grouping states and their evolution. 
The grouping phenomenon has also been revealed in a real-world example of the subsystem of 27 stocks in 
the Shanghai Stock Market's Steel Plate. Our work demonstrates that complex systems following the MG 
rules can spontaneously self-organize themselves into certain divided states, and our model represents a basic 
mathematical framework to address this kind of phenomena in social, economical, and even political systems. 

PACS numbers: 02.50.Le, 89.75.Hc, 87.23.Ge 



I. INTRODUCTION 

The Minority Game (MG) was originated from the El Parol 
bar problem in game theory first conceived by Arthur in 1994 
111], where a finite population of people try to decide, at the 
same time, whether to go to the bar on a particular night. Since 
the capacity of the bar is limited, it can only accommodate a 
small fraction of all who are interested. If many people choose 
to go to the bar, it will be crowded, depriving the people of the 
fun and thereby defying the purpose of going to the bar. In this 
case, those who choose to stay home are the winners. How- 
ever, if many people decide to stay at home then the bar will be 
empty, so those who choose to go to the bar will have fun and 
they are the winners. Apparently, no matter what method each 
person uses to make a decision, the option taken by majority 
of people is guaranteed to fail and the winners are those that 
choose the minority one. Indeed, it can be proved that, for the 
El Parol bar problem there are mixed strategies and a Nash- 
equilibrium solution does exist, in which the option taken by 
minority wins |2]. A variant of the problem was subsequently 
proposed by Challet and Zhang, named as an MG problem 
yfl, where a player among an odd number of players chooses 
one of the two options at each time step. Subsequently, the 
model was studied in a series of works |4-22]. In physics, 
MG has received a great deal of attention from the statistical- 
mechanics community, especially in terms of problems asso- 
ciated with non-equilibrium phase transitions ll23l - l25ll . 

In the current literature, the setting of MG is that there is 
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a single resource with players' two possible options (e.g., in 
the El Parol bar problem there is a single bar and the options 
of agents are either going to the bar or not), and an agent 
is assumed to react to available global information about the 
history of the system by taking on an alternative option that 
is different than the current one it is taking. The outstand- 
ing question remains of the nonlinear dynamics of MG with 
multiple resources. The puipose of this paper is to present a 
class of multi-resource MG models. In particular, we assume 
a complex system with multiple resources and, at any time, 
an individual agent has k > 1 resources/strategies to choose 
from. We introduce a parameter p, which is the probability 
that each agent responds based on its available local infor- 
mation by selecting a less crowded resource in an attempt to 
gain higher payoff. We call p the minority-preference proba- 
bility. We find, strikingly, as p is increased, the phenomenon 
of grouping emerges, where the resources can be distinctly 
divided into two groups according to the number of their at- 
tendees. In addition, the number of stable pairs of groups also 
increases. We shall show that the grouping phenomenon plays 
a fundamental role in shaping the fluctuations of the system. 
The phenomenon will be demonstrated numeiically and ex- 
plained by a comprehensive analytic theory. An application 
to the analysis of empirical data from a real financial market 
will also be illustrated, where grouping of stocks (resources) 
appears. Our model is not only directly relevant to nonlinear 
and complex dynamical systems, but also applicable to social 
and economical systems. 

Our multi-resource MG model is presented in Sec. jllj The 
emergence of grouping phenomenon is demonstrated in Sec. 
nil] A general theory is developed in Sec. |IV] to elucidate 
the dynamics of the emergence and evolution of the strategy 
groups. Concluding remarks are presented in Sec. jVT] 
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II. MULTI-RESOURCE MINORITY GAME MODEL 

We consider a complex, evolutionary-game type of dynam- 
ical system of N interacting agents competing for multiple 
resources. Each agent will chose one resource in each round 
of the game. And, each resource has a limited capacity, i.e., 
the number of agents it can accommodate has an upper bound 
ric- There are thus multiple strategies {s — 1, 2, ■ ■ ■ , k, where 
k is the maximum number of resources/strategies) available to 
each agent. On average, each strategy can accommodate N/ k 
agents, and we consider the simple case of ric = N/k. Let 
Ug be the number of agents selecting a particular strategy s. If 
Us < Tic, the corresponding agents win the game and, conse- 
quently, s is the minority strategy. However, if Ug > ric, the 
associated resource is too crowded so that the strategy fails 
and the agents taking it lose the game, which defines the "ma- 
jority strategy." The optimal solution to the game dynamics is 

thus Us — He- 

In a real-world system, it is often difficult or even impossi- 
ble for each agent to gain global information about the dynam- 
ical state of the whole system. It is therefore useful to intro- 
duce the concept of local information network in our multiple- 
resource MG model. At each time step, with probability p, 
namely the minority-preference probability, each agent acts 
based on local information that it gains by selecting one of 
the k available strategies. In contrast, with probability 1 — p, 
an agent acts without the guidance of any local information. 
For the minority-preference case, agent i has d neighbors in 
the networked system. The required information for i to react 
consists of all its neighbors' strategies and, among them, the 
winners of the game, i.e., those neighboring agents choosing 
the minority strategies at the last time step. Let 11 = {s„i} be 
the set of minority strategies for i's winning neighbors, where 
a strategy may appear a number of times, if it has been chosen 
by different winning neighbors. With probability p, agent i 
will chose one strategy randomly from 11. Thus, the probabil- 
ity Ps for strategy s to be selected is proportional to the times 
it appears in 11, i.e., Pg = Ns/Card{ll), where Card{ll) is 
the number of elements in 11 and Ng the times strategy s ap- 
pears in n. If n is empty, i will randomly select one from 
the k available strategies. While, for the case that an agent se- 
lects a strategy without the guidance of any local information 
with probability 1 — p, it will either choose a different strategy 
randomly from the k available ones with mutation probabil- 
ity m 1 10], or inherit its strategy from the last time step with 
probability 1 — to. 



III. NUMERICAL RESULTS 

As a concrete example to illustrate the strategy-grouping 
phenomenon, we set A; = 5. Figures [Ua-c) show time se- 
ries of Ug, the number of agents selecting each strategy s, for 
p = 0, 0.45, and LO, respectively. For Fig. [Tfa) where p = 0, 
an agent makes no informed decision in that it changes strat- 
egy randomly with probability m but stays with the original 
strategy with probability 1 — m. In this case, n^'s appear ran- 
dom. For the opposite extreme case of p = 1 [Fig. [TJc)], each 



agent makes well informed decisions based on available lo- 
cal information about the strategies used by its neighbors. In 
this case, the time series are quasiperiodic (a detailed analy- 
sis will be provided in Sec. llVb . For the intermediate case of 
p — 0.45 [Fig. [Tfb)], agents' decisions are partially informed. 
In this case, an examination of the time series points to the oc- 
currence of an interesting grouping behavior: the 5 strategies, 
in terms of their selection by the agents, are divided into two 
distinct groups gi and 52 that contain kg^ ~ 2 and kg.^ = 3 
resources, respectively. The time series associated with the 
smaller group exhibit larger fluctuations about its equilibrium. 

To better characterize the fluctuating behaviors in the time 
series Ug, we calculate the variance al = ([«s(0 — N /k]"^) as 
a function of the system parameter p, where (■) is the expec- 
tation values averaged over a long time interval, as shown in 
Fig. |2]on a logarithmic scale. We observe a generally increas- 
ing behavior in a1 with p and, strikingly, a bifurcation-like 
phenomenon. In particular, for p < pf,, where p^ is the bifur- 
cation point, (j^s for all strategies assume approximately the 
same value. However, for p > pb, there are two distinct val- 
ues for a^, signifying the aforementioned grouping behavior 
[Fig- nib)]. From Fig. |2] we also see that, after the bifurca- 
tion, the two branches of ct^ are linear (on a logarithmic scale) 
and have approximately the same slope a, suggesting the fol- 
lowing power-law relation: (7^. — bg^p"", for i = 1,2, where 
log ) and log {hg^ ) are the intercepts of the two lines in 
Fig. |2] We thus obtain 

4^ = ^. (1) 

In Sec. IIVI we will develop a theory to explain the relations 
among the variances of the grouped strategies and to provide 
formulas for the amplitudes of the time series in Fig. [Tjand 
the sizes of the groups (denoted by kg^ and kg^, respectively). 
Specifically, our theory predicts the following ratio between 
the variances of the two bifurcated branches: 

^ = (^)^ (2) 

which is identical to the numerically observed ratio in Eq. ([T]i, 
with the additional prediction that the strategies in the group 
of smaller size exhibit stronger fluctuations since the corre- 
sponding value of CT^ is larger Overall, the emergence of the 
grouping behavior in multiple-resource MGs, as exemplified 
in Fig. 12] resembles a period-doubling like bifurcation. While 
period-doubling bifurcations are extremely common in non- 
linear dynamical systems, to our knowledge, in complex game 
systems a clear signature of such a bifurcation had not been 
reported previously. 

A careful examination of the time-series for various p 
has revealed that the strategy-grouping processes has already 
taken place prior to the bifurcation point pb in the variance 
(j"^, but all resulted grouping states are unstable. Take as an 
example the 5-strategy system in Fig. [T] In principle, there 
can be two types of pairing groups: (1, 4) and (2, 3). For any 
grouping state, the following constraint applies: 

^91 ^ ^92 ~ ^91 ■ ^92 7^ (■^) 
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FIG. 1 : (Color online.) (a-b) For a square-lattice system of A'^ = 10000 agents, time series of the number of agents selecting each of the k — 5 
available strategies for p — 0, 0.45, and 1 .0, respectively. The probability of random alteration of strategy in the case of complete lack of local 
information is set to be m = 1. A strategy grouping behavior can be seen in (b), where the whole strategy set is broken into two groups: one 
of two and another of three strategies. 
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FIG. 2: For the same system in Fig. [T] variance of the time series 
Ua as a function of p on a logarithmic scale. A period-doubling like 
bifurcation occurs, at which a grouping behavior emerges. Inset is 
the same plot but on a linear scale. 



There are in total fc/ 2 (if A: is even) or (fe — l)/2 (if fc is odd) 
possible grouping states for the system with k available strate- 
gies. However, the grouping states are not stable for p < pi,. 
What happens is that a strategy can remain in one group but 
only for a finite amount of time before switching to a differ- 
ent group. Assume that the sizes of the original two pairing 
groups are (fc^^ and kg^), respectively. The sizes of the new 
pair of groups are thus {kg^ ± 1 and kg^ T 1)> as stipulated 
by Eq. Associated with switching to a different pair of 
groups, the amplitudes of the time series Us for each strategy 
also change. As the bifurcation parameter p is increased, the 
stabilities of different pairs of grouping states also change. At 
the bifurcation point pi,, one particular pair of groups becomes 
stable, such as the grouping state (2,3) in Fig. |2] 

The bifurcation-like phenomenon and the emergence of 



various strategy-grouping states are general for multiple - 
resource MG game dynamics. For example. Fig. |3]shows 
as a function of p for a system with fc = 16 available strate- 
gies. There are in total 8 possible grouping states, ranging 
from (8,8) to (1, 15). As p is increased, the grouping states 
(8, 8), (7, 9), (6, 10), (5,11) and (4, 12) become stable one af- 
ter another, as can be seen from the appearance of their corre- 
sponding branches in Fig. [3] The behavior can be understood 
theoretically through a stability analysis (Sec. lIVI i. 

Another phenomenon revealed by Fig. |3]is the merging of 
bifurcated branches. For example, as p is increased through 
about 0.8, the grouping states disappear one after another in 
the reverse order as they initially appeared. This can also be 
understood through the stability analysis (Sec. llVt . 



IV. THEORY 

Here we develop an analytic theory to understand the emer- 
gence, characteristics, and evolutions of the strategy groups. 



A. Relationship among variance, amplitude and group size 

In general, for a multiple-resource MG system of N agents, 
as the parameter p is increased so that agents become more 
likely to make informed decision for strategy selection, the 
available strategies can be divided into pairs of groups. The 
example in Fig. [Ttb) presents a case where there are two 
distinct strategy groups gi and g2, which contain fcg^ and 
fcg2 Strategies, respectively, where fcg^ + kg^ — k. For Fig. 
[Tfb), we have fcg^ < kg^. The strategies belonging to the 
same group are selected by approximately the same number 
of agents, i.e., the time series ns{t) for strategies in the same 
group are nearly identical. During the time evolution, a strat- 
egy s can switch iteratively from being a minority strategy 
[ns{t) < ric = N/k] to being amajorityone [ns{t+l) > ric]- 
In particular, as shown in the schematic map in Fig. |4] for the 
strategy in group gi denoted by Si, and the strategy in group 



FIG. 3: (Color online.) For a multiple-resource MG system of A'^ = 10000 agents on a square lattice and fc = 16 available strategies, (a) 
numerically obtained bifurcation-like behavior that leads to the emergence of various pairs of grouping states [e.g, (8,8), (4,12), etc.], (b) 
predicted bifurcation from mean-field theory (Sec. IIV) . and (c) results from an improved mean-field theory (Sec. II V) . The probability of 
random selection in the absence of local information is set to be m = 1. 
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FIG. 4: For the case of two groups [e.g.. Fig. [TJb)], schematic illus- 
tration of time series Us^ in group gi (blue circles) and Us^ in group 
(72 (red triangles), where various quantities such as the amplitude A 
are labeled. 



is closed so that the number of agents is a constant, we have 

fcffi • (t) + kg, ■ (t) = N. (5) 
Thus, for two consecutive time steps, we have, 

kg, ■ {{us^) - As,) + kg, ■ {{n,^) + As,) = N. (6) 
Substituting Eq. (|4]i into Eq. (|6]), we have 



,N 



,N 



kg, ■ (— + An,, +AsJ + kg, ■ (— - Aus^ - As 

N N 
kg, ■ (— + Aus^ - + kg, ■ (— - Aus^ +As;) = N, 

from which we obtain the relations between As- and As,, 
An,, and An, . as 



k 



92 



^92 



kg, AtLs, kg. 



(7) 



f/2 denoted by Sj, if ns-{t) < N/k < ns,{t), we will have 
nj. (i + 1) > N/k > Us, {t + 1). In addition, the time seiies 
Us reveals that the average numbers of agents for strategies 
Si and Sj, denoted by (n,. } and (n,^ ) (the blue dash line and 
red dot line in Fig. |4]i, respectively, are not equal to N/k (the 
black solid line in Fig. |4|i. In fact, we have 



At 



At 



(4) 



with {ns/) > N/k, and (n,^.) < N/k. Here, An,^ is the 
absolute value of the difference between (us^) and N/k. The 
number of agents n^^ {t) [or Us, (<)] typically fluctuates about 
the equilibrium (n,;) [or (n,^. )] with amplitude As- [or As,], 
as shown in the schematic map in Fig. |4] 

Based on the numerical observations, we can argue that 
the strategy grouping phenomenon is intimately related to the 
fluctuations in the time series ris. Assuming the MG system 



We see that the fluctuations of the time series are closely re- 
lated the grouping of the strategies. 

From the definition of cr^, we can write the variances of n^^ 
and n, . as 



{As^ + Aus^)"^ + [As^ -Aus^)^ 
2 

{As^ ~Ans^)^ + {As^ +AnsS' 



Ai.+Ani 



Al + Anl 



Using Eq. (|7]i, we obtain 



^91 



(8) 



As shown in Fig. |2] the ratio of the variances of group gi and 
92 from the simulation agree very well with Eq. (|8). 
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B. Mean-field theory 

We develop a mean-field theory to understand the fluctua- 
tion patterns of the system. To be concrete, we still treat the 
case of two distinct groups. Consider strategy Si that belongs 
to group gi and assume that Si is the majority strategy at time 
t = to, i.e., rti^' > N/k. According to the mean-field approx- 
imation, at the next time step t 
nip choosing strategy Si is 



to + 1, the number of agents 



The process Ug 



.(1) 



,(2) 



thus occurs iteratively. From 



Eqs. ( fTTT i and ( fT3] l, we can get the number of agents for one 
given strategy at any time t. In particular, denoting a = (1 — 

p){l — m), (3 = N{1 —p)m/k, and = Np/kg^, we obtain 
the iterative dynamics for agents selecting strategy Si as 



(14) 



.(1) 



,(0) 



,(0) 



(wi + W2) + W3 



(0)/ 



(iV~n(°))(l-p)mi 



(9) 



where 



Wi — nf.^p, 
,(0) 



(fc-1) 



■013 = {N - n(°))(l -p)m 



k' 



(10) 



Here, wi and W2 together are the number of agents abandon- 
ing strategy Si (or the flow out of strategy Si). That is, of the 
nf.^ agents, wi agents will act based on local information by 
selecting a minority strategy different than Si for the next time 
step to + 1. At the same time, there will be W2 agents acting 
without local information by choosing randomly one of the 
other fc — 1 strategies. The quantity W3 represents the flow 
into Si from the remaining TV — nif' agents. These agents will 
mutate randomly to switch their strategies to Si without any 
local information. We thus have 



,(1) 



,(0) 



';>{l-p){l-m)+N{l-p)m 



(11) 



Suppose ni; < N/ k. Then Si and the other kg-^ — 1 strategies 
in gi are the minority strategy, and the agents selecting those 



strategies win the game att = to 
time to + 2 can be written as 



1. The time series n. 



(2) 



at 



.(2) 



.«-nW[p%-^ + (l-p)mi^] 



+ {N ~ nmp— + {1 - p)m-] 



(W4 + W5) + {W6 + W'j), 



(12) 



where and W5 stand for the flows out of strategy s;, while 
wg and w-j represent the flows into Sj from those N — n'fj 
agents on other strategies. We have 



7(2) ^„(i)n - 



1 



1 



ni\>{l-p){l - to) + N{1 -p)m- + Np—, (13) 

k kg^ 

where tis/ is larger than N/k, Si and all other strategies in 
group gi will be the majority strategy again, as at time to- 



,(2a 



,(2a+l) 



an 



(2a) 



f3, 



where a G R and ni\ stands for the number of agents at 
t = to +t'. Carrying out the iterative process in Eq. (Il4l . we 
obtain 

=a^nfj +al3 + /3 + jg,, 
nif = a^nW + a/S + ajg^ + {3. 

For the case where ns{t) exhibits stable oscillations, i.e., the 
system is in a stationary state, we have. 



(2a) 



nW = ...=ni^«+i). 

We thus obtain the values of ni^"' and ni^"^^' as a function 
of the probability p, mutation probability to, and the grouping 
parameter fc^^, and fc: 

(2a) _ al3 + l3 + ig. 

Si ~ 1 2 ' 

„(2a+l) ^ «/3 + QTgi + /3 



1 



(15) 



From Tis, we can get the expression of the ampUtude of the 
fluctuation, the mean value (n^ . ), and its difference from N/k 
as. 



2(1 + a) ^ 



An. 



IK)-7v/fc| = | 



2(1 -a)' 

fc7g, -2Np _ Np\k/kg, 

2k{l-a) ' 



2fc(l 



In the above derivation, we have assumed that Si belongs to 
group gi, and obtained expressions of the ns-. As-, (n-s.), and 
Aus- ■ Similarly, using Eq. (|5]l, we can calculate the time se- 
ries nsj , the number of agents selecting the strategy sj belong- 
ing to group g2, and the corresponding characterizing quanti- 
ties. 

Alternatively, following the steps similar to those from Eqs. 
(|9]l to ( fT4l i. we can write the recurrence formula for rig as 



nl^-^) = an 



nif+^)=an(^'^)+/3 + 7,., 
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FIG. 5: (Color online.) Comparison of the time series of rij,. in 
group gi (and Usj in group 32) from the mean-field formula [Eq. 
J15l >1 and modified mean-field formula [Eq. dlSH. and from simula- 
tions. The upper and lower figures show the results for the grouping 
state (8,8) and (7,9), respectively. The left panels are the results 
from mean-field theory (MF) and from simulations on fully con- 
nected networks (FCN). The right panels are the results from our 
modified mean-field theory (MMF) and simulations on the square lat- 
tice (SL). The system parameters are A*' = 10000, k — 16, p — 0.45, 
and m — 1. 



where jg^ — Np/kg^. For the case of stable strategy, we get 
the corresponding number of agents as 



(2a) ^ a/3 + 0792 + P 
1 - a2 

(2a+l) _ + + 7g2 

" l-a2 ■ 



(16) 



We find that the values of Us obtained from Eq. (fT6] l agree 
with those from Eq. (|5]l very well (derivation and data not 
shown). In addition, the expressions of Ag^, {ngj), and Aug^ 
are identical to those associated with Sj, with the quantity 7^^ 
replaced by jg^ . 

The mean-field theory is ideally suited for fully connected 
networks. Indeed, results from the theory and direct simu- 
lations agree with each other very well, as shown in Fig. |5] 
However, in real-world situations, a fully connected topol- 
ogy cannot be expected, and the mean-field treatment wiU no 
longer be accurate. For example, we have carried out sim- 
ulations on square-lattice systems and found noticeable de- 
viations from the mean-field prediction. To remedy this defi- 
ciency, we develop a modified mean-field analysis for MG dy- 
namics on sparsely homogeneous networks (e.g., square lat- 
tices or random networks). 



C. Mean-field theory for MG dynamics on sparsely 
homogeneous networks 

Due to the limited number of links in a typical large-scale 
network, it is possible for a failed agent to be surrounded by 



agents from the same group (who will likewise fail the game). 
In this case, the failed agent has no minority strategy to imitate 
(set n is empty) and thus will randomly select one strategy 
from the k available strategies. Taking this effect into account, 
we can modify the mean-field approximation in Eqs. (|9]l and 
( fT2l i as 



,(2) 



(0) 



{wi + W2) + W3 + Ws, 



(17) 



with the two modified terms given by 

Wq^{N - kg,n[]^)pr]g^{j - -5-), 

where rjg^ is the probability for one agent in group to be 
suiTounded by agents from the same group. The quantity 
stands for the flow from the failed agents in group gi [the 
number is Ng^ = fc^^ni"''], who react to the information [the 
number is Ng^p} but with no winner surrounded to supply the 
optional minority strategy [the number is Ng-^prjg], and thus 
select Si with probability 1/fc. The quantity wg represent two 



factors: (1) the rii = {N — kg 



)P'ir~''^92 failed agents in 



group g2 who should have flow into Si [i.e. wq] but are held 
back because they are surrounded by agents in g2, and (2) 
the failed agents in group g2 who are surrounded by agents 
in g2 and thus select Si with 1/fc probability, the number of 
which is n2 = {N — kg^rLi})prig^^. Apparently, we have 

From Eq. ( fTTj l. we obtain 



(2a) ^ + P + 

(2a+l) ^ + +/^ 

1 - a' a" 



(18) 



As. 
An.,, = 



2(1 + 

2) + 73,(1 + 
2{l -a'g^a';,^) 

| /3«+<)+7;,(l + 0+2/3 N 
' 2{l-a'gy^J fc' 



where the parameters are 

< ^ (l-J3)(l-"^) + 



"si (1 - P)(l - n^) + kgiPVg2 ( 
7;=^M^+%.(^-^)]- 



The equation set dTsl l represents the modified mean-field de- 
scription of the time series associated with the stable strate- 
gies in the game system supported on sparsely homogeneous 
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FIG. 6: (Color online.) Phase diagram indicating the stability region 
of strategy-grouping states. The stable parameter region of one given 
grouping state with fcg^ — k' is at the upper right of the correspond- 
ing curve of k' . 



networks. The density of agents in is denoted by pg^ = 
Ng^/N. For the case where agents from different groups 
are well mixed in the network, the probability for one given 
agent in to meet with agents in gy is Pg^g^ — Pg^ (for 
x^y — 1,2). If the average degree of the network is d, the 
probability that one agent from g^ is surrounded by agents 
from the same group is rjg^ = [Pg^g^Y'' = Pg^- From the 
simulation on the square lattice system where each agent has 
d = 4 neighbors, we observe a quite weak effect of clustering 
of agents in gi or g2, so rjg^ sa Pg^ — (Ng^/N)'^. Based on 
the quantity rjg^ , analytic prediction of the time series in the 
modified mean-field (MMF) theory can be obtained, as shown 
in Fig. |5]for a square-lattice system. We observe a good agree- 
ment with simulation results. Simulations on homogeneous 
networks of different values d have also been carried out, with 
results in good agreement with the prediction from the modi- 
fied mean-field theory. 



Using Eqs. ( fTsT i and ( fTSb , we get 

— m — p + pm) 
2 — m — p + pm 
k 



kg^ > 



kq2 ^ "T — 

2 — m — p + mp 



(19) 



where and ^2 are continuous functions of the parameters p 
and 771, and ^1+^2 = k. The two inequalities in Eq. (l% 
are nevertheless equivalent to each other Figure |6]presents 
a phase diagram in the parameter space, where the curves of 
^1 {p, m) = k' for fc' = 8, 7, • • • , 1 are shown. The necessary 
condition for the strategy-grouping state with kg^ = k' to be 
stable is that the the parameters p and 777, are in the upper- 
right region of the curve 1^1(^,777) = k' . For certain value 
of 777,, only when p > pb{rn,k') will the state of kg^ = k' 
be stable. While the value of pb from simulation is different 
from the theoretical value pi,{m,k'), our mean-field theory 
does provide a qualitative explanation for the phenomenon in 
Fig. ISa), where more branches of smaller strategy-grouping 
states become stable as p is increased. 

We have also seen in Fig. |2] and Fig. [Sja) that, as p ap- 
proaches 1, the bifurcated branches of different grouping state 
merge together For the case of even k, Ug fluctuates stably in 
the grouping state with fcg^ — kg^ = fc/2. While, for the case 
of odd k. Us fluctuates quasiperiodicly [see Fig. [Uc)]- Actu- 
ally, the grouping state always switches between {kg + 1, kg) 
and (fcg, fcg + 1), with kg ~ (fc— l)/2. We can understand the 
instability and merging of grouping states from Fig. [TJc), and 
the schematic map in Fig. |4] as follows. As p is increased to 1, 
Aus- and Aus- increase and become comparable to the am- 
plitudes and , respectively. Namely, the attendances 
r7s of strategies can be very close to ric — N/k . In case 
that Ug of one strategy does not get cross Uc because of noise, 
i.e., it acts as minority (or majority) strategy twice, then, the 
fluctuation of Ug, as well as the grouping state is changed. 



V. A REAL- WORLD EXAMPLE: EMERGENCE OF 
GROUPING STATES IN FINANCIAL MARKET 



D. Stability of strategy grouping states 

Our mean-field treatment yields formula characterizing 
the stable oscillations associated with the grouping state 
(fcgjjfcg^), which include the variance ratio of the pairing 
groups [Eq. (O] and time series 77^ [Eqs. ( fTSl ) and (fTSlll. How- 
ever, from the simulation result shown in Fig. |3] we see that 
not all the grouping states are stable in the parameter space. 
As p is increased, the strategy-grouping state of the smaller 
fcgj becomes stable, and the coiTesponding branch appears. 
It is therefore useful to analyze the stability of the grouping 
state. 

In our treatment we have assumed ni^^ > N/ fc and 
N/k. Then, the necessary condition for the grouping state 
(fcgj, fcgj) to become stable is 11)^} < N/k and 77!^'' > N/k. 



The financial market is a representative multi-resource 
complex system, in which many stocks are available for in- 
vestment. We analyze the fluctuation of the stock price from 
the empirical data of 27 stocks in the Shanghai Stock Mar- 
ket's Steel Plate between 2007 and 2010. We regard the 27 
stocks, which are derived from the iron and steel industry, 
as constituting a MG system with k — 27 resources, where 
the agents selecting the resources correspond to the capitals 
invested. This system is open in the sense that capital typi- 
cally flows in and out, which is the main difference from our 
closed-system model. In particular, given the time series Xi {t) 
of the daily closing price of stock i, the daily log-return is 
defined as Ri{t) = In Xi{t) — In Xi{t — 1). The average re- 
turn of the 27 stocks at time t, denoted by {Ri{t)), signifies a 
global trend of the system at t, which is caused by the change 
in the total mount of the capital in this open, 27-stock sys- 
tem. However, when we analyze the detrended log-returns 
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FIG. 7: (Color online.) Correlation matrix of log-returns for 27 stocks in Shanghai Stock Market's Steel Plate, (a) Original matrix C, (b) 
matrix C' ordered according to the eigenvector for the maximum eigenvalue of C, and (c) reordered matrix C" from C' with respect to d 
within each group. The Ci values of the stocks are shown for each matrix. 



R[{t) — Ri{t) — {Ri{t)), the system resembles a closed sys- 
tem, as our model MG system. We shall demonstrate that the 
strategy-grouping phenomenon occurs in this real-world sys- 
tem. 

We calculate the Pearson parameter of each pair of the 
detrended log -returns R[ [t) and R'^ (t), which leads to a fc x fc 
correlation matrix C, as shown in Fig. Ha). In terms of the 
eigenvector associated with the maximum eigenvalue of ma- 
trix C, we rank the order of the stocks and obtain the matrix 
C', as shown in Fig. |3b). The striking behavior is that the ma- 
trix is apparently divided into 4 blocks, a manifestation of the 
grouping phenomenon. In particular, the matrix elements Cij 
among the first 15 stocks and those among the remaining 12 
stocks are generally positive, but the cross elements between 
stocks in the two groups are negative. It is thus quite natural 
to classify the first 15 stocks as belonging to group gi and the 
remaining 12 to group 52- We can then write the matrix in a 
block form as 

Q! = ( ^9291 ^3232 j 

where the elements of Cg^g^ and Cg^g^ are positive, and those 
of Cgjgj and Cg^g-^ ciTS uegatlvc. The phenomenon is that 
the 27-stock system has self-organized itself into a (12, 15) 
grouping state, which is a natural consequence of the MG dy- 
namics in multi-resource complex systems. 

For one given stock i, the mean absolute correlation is 

fc 

Ci = ^\cij\/{k-l), j^i. 
j=i 

This parameter reflects the weight of the stock in the system. 
If Ci — > 0, oscillations of stock i are contained in the noise 



floor. In this case, there is no indication as to whether this 
stock belongs to group gi or 52- The larger the value of 
the less ambiguous that the stock belongs to either one of the 
two groups. From the value of q ranked in the same order 
as in C', we can see that the boundary of the two groups is 
the stock with minimum Ci . Thus Ci can be considered as the 
characteristic number to distinguish different groups. We can 
also reorder the matiix C' according to q within group gi and 
g2, respectively. This leads to the matrix C", as shown in Fig. 
|3c), further demonstrating the grouping phenomenon. 



VI. CONCLUSIONS AND DISCUSSIONS 

Minority game, since its invention about two decades ago, 
has become a paradigm to study the social and economical 
phenomena where a large number of agents attempt to make 
simultaneous decision by choosing one of the available op- 
tions [2J. In the most commonly studied case of a single avail- 
able resource with players' two possible options, agents tak- 
ing the minority option are the guaranteed winners. Various 
minority game dynamics have also received attention from the 
physics community due to their high relevance to a number 
of phenomena in statistical physics. It has become more and 
more common in the modern world that multiple resources 
are available for various social and economical systems. If 
the rule still holds that the winning options are minority ones, 
the questions that naturally arise are what type of collective 
behaviors can emerge and how they would evolve in the un- 
derlying complex system. Our present work aims to address 
these questions computationally and analytically. 

The main contribution and findings of this paper are the 
following. Firstly, we have constructed a class of spatially 
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FIG. 8: (Color online.) (a) The analytical time series ns(t) of period-tripling in the form of 2 groups, and (b) the corresponding mapping 
from Hs - (t) to ris . (t + 1) for the strategy Si in group gi [see Eq. (IIVH. The system is of fc = 16, p — 0.25, m — 0.2, and the grouping 
state (fegj , ) is (3, 13), respectively, (c) The time series Us (t) of the stable grouping state with period-tripling in the form of 3 groups from 
simulation on FCN. The system is of fc = 16, p = 0.88, m — 1, and the grouping state (fcg^ , kg^ , kg^) is (5, 5, 6). 



extended systems in which any agent interacts with a finite 
but fixed number of neighbors and can choose either to follow 
the minority strategy based on information about the neigh- 
boring states or to select one randomly from a set of avail- 
able strategies. The probability to follow the local minority 
strategy, or the probability of minority preference, is a key 
parameter determining the dynamics of the underlying com- 
plex system. Secondly, we have carried out extensive nu- 
merical simulations and discovered the emergence of a strik- 
ing collective behavior: as the minority-preference probabil- 
ity is increased through a critical value, the set of available 
strategies/resources spontaneously break into pairs of groups, 
where the strategies in the same group are associated with a 
specific fluctuating behavior of attendance. This phenomenon 
of strategy-grouping is completely self-organized, which we 
conjecture is the hallmark of MG dynamics with multiple 
resources. Thirdly, we have developed a mean-field theory 
to explain and predict the emergence and evolution of the 
strategy-grouping states, with good agreement with the nu- 
merics. Fourthly, we have examined a real-world system of 
a relatively small-scale stock-trading system, and found un- 
equivocal evidence of the grouping phenomenon. Our results 
suggest grouping of resources as a fundamental type of collec- 
tive dynamics in multiple-resource MG systems. Other real- 



world systems for which our model is applicable include, e.g., 
hedge-fund portfolios in financial systems, routing issues in 
computer networks and urban traffic systems. We expect our 
model and findings to be not only relevant with statistical- 
physical systems, but also important to a host of social, eco- 
nomical, and political systems. 

Additionally, the double-grouping (or paired grouping) 
period-doubling like bifurcation are not the exclusive mode 
for the grouping of resources. We have also observed 
the period— 3 double-grouping [see Fig. [8ja)(b)], and the 
period— 3 triplet-grouping bifurcation phenomena [see Fig. 
|8jc)]. Further research about these phenomena are of great 
interest and is valuable for the understanding of chaotic and 
bifurcation in social systems. 
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